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Abstract 

We consider a statistical system in a planar wedge, for values of the bulk parameters cor¬ 
responding to a first order phase transition and with boundary conditions inducing phase 
separation. Our previous exact field theoretical solution for the case of a single interface is 
extended to a class of systems, including the Blume-Capel model as the simplest represen¬ 
tative, allowing for the appearance of an intermediate layer of a third phase. We show that 
the interfaces separating the different phases behave as trajectories of vicious walkers, and 
determine their passage probabilities. We also show how the theory leads to a remarkable 
form of wedge covariance, i.e. a relation between properties in the wedge and in the half 
plane, which involves the appearance of self-Fourier functions. 


Melfino@sissa.it, Mlessio.squarcini@sissa.it 


1 Introduction 


Fluid interfacial phenomena at boundaries form an important chapter of statistical physics 
and are studied experimentally, theoretically and numerically (see p [21 [31 n 0 El [3 0 [g do] 
for reviews). On the theoretical side, the exact results obtained for the lattice Ising model 
in two dimensions muniE] provided an important benchmark for approximated or heuristic 
approaches, but proved too difficult to extend to other universality classes. Only recently it has 
been shown that phase separation and the interfacial region in planar systems can be described 
exactly for the different universality classes [lainiiis] relying on low energy properties of two- 
dimensional field theory m- This new approach also allowed the exact solution m of the 
longstanding problem of phase separation in a wedge, which received much attention [lEiiiniEni 
[211 [22l [23l [231 [231 [26] as the basic example of the effect of the geometry of the substrate on the 
adsorption properties of a fluid. Particular interest was attracted by emergent relations between 
adsorption properties in the wedge and those on a flat substrate. Observed at the macroscopic 
level [IB] and successively referred to as properties of “wedge covariance”, this type of relations 
resisted a derivation within a statistical mechanical framework. It was one of the results of m 
to show for the planar case how wedge covariance follows from the relativistic invariance of the 
quantum field theory associated to the universality class, which in turn reflects the homogeneous 
and isotropic nature of the fluid. 

The wedge problem has been considered so far for the case in which the universality class and 
the boundary conditions lead to the separation of two phases o and b. Here we develop our exact 
field theoretical approach to study the case in which a macroscopic bubble of a third phase c 
forms in between the two favored by the boundary conditions. Establishing whether a third phase 
will intrude between phases a and b forming a macroscopic intermediate layer or just microscopic 
droplets at the interface is a main question of wetting physics. It was shown in |13l [T3] that 
in two dimensions the answer is determined by the spectrum of elementary excitations of the 
underlying field theory, which is known for the different universality classes. A model which 
leads to the macroscopic (wetting) layer of the third phase and we will consider in this paper is 
the g-state Potts ferromagnet m at its first order transition point, at which the q ferromagnetic 
phases (two of which correspond to phases a and b) coexist with the disordered phase (which will 
play the role of phase c). In two dimensions the first order transition corresponds to T = Tc for 
(7 > 4 |28j : strictly speaking, our field theoretical description is exact for the scaling limit —)■ 4+ 
|29j . but is expected to remain quantitatively meaningful up to g ~ 10, where the correlation 
length is still much larger than lattice spacing. Our derivation applies also to g < 4 provided we 
allow for the possibility of vacant sites (dilution); then coexistence of the disordered with the 
ordered phases is recovered above a critical value of dilution for a value Tc of the temperature. 
For g = 2 one obtains a dilute Ising model, also known as Blume-Capel model, for which the 
wetting character of the disordered phase has been investigated numerically [301 ED Eg. 

We will work for values of the boundary parameters such that the inner phase is not adsorbed 
on the boundary. This means that the third phase is separated from the other two by two 


1 


interfaces which fluctuate between the two boundary conditions changing points. We will see 
how this picture emerges within the field theoretical framework and will determine the passage 
probabilities for the interfaces, finding in particular that they randomly fluctuate with the 
constraint of avoiding each other and the boundary, i.e. that they correspond to trajectories of 
so-called “vicious” walkers [33]; in this way we determine the passage probabilities of vicious 
walkers in a wedge. Concerning the issue of wedge covariance, it turns out to acquire additional 
interest in the case of two interfaces, with a surprising interplay between physical considerations 
in momemtum space and mathematical realization of the condition of impenetrability of the 
wedge. 

The derivations are exact and apply to the case of a shallow wedge, that for which the 
results are universal, in the sense that they do not depend on the specific values of the boundary 
parameters, as long these are in the range which does not bind the interfaces to the boundary. 

The paper is organized as follows. In the next section we recall the setting and the results 
of [T7|. This will put us in the condition of developing the theory for the case of two interfaces 
in section 3. The final section is then devoted to summary and comments. 

2 Two phases in a wedge 

We start with the characterization of the statistical system in absence of boundaries, i.e. on the 
infinite plane. The system is considered at a first order phase transition point, where different 
phases, that we label by an index a = 1, 2,... , n, have the same free energy and can coexist at 
equilibrium. At the same time the system is supposed to be close to a second order transition 
poin10, in such a way that the correlation length is much larger than microscopic scales and 
a continuous description is allowed. For homogeneous and isotropic systems this continuous 
description is provided by a Euclidean field theory with coordinates (x, y) identifying a point on 
the plane. This field theory in turn corresponds to the continuation to imaginary time t = iy of 
a quantum field theory in one space dimension with coordinate x. The degenerate phases of the 
statistical system are in one-to-one correspondence with degenerate vacua |0a) of the associated 
quantum theory. We denote by a{x,y) the order parameter field, and by {a)a = (0a|<7(x,y)|0a) 
the value of the order parameter in phase a. For a generic field we have 

^>(x, y) = e^^-“^^>(0,0)e-^^+“^ , (1) 

with the Hamiltonian H and momentum operator P of the quantum system acting as generators 
of time and space translations, respectively; the vacuum states carry zero energy and momentum. 

As usual in presence of degenerate vacua in (H-l) dimensions (see (Ej), the elementary 
excitations correspond to kinks \Kab{0)) which interpolate between two different vacua |0a) and 
|0fe), and whose energy and momentum satisfy the relativistic dispersion relation 

(e,p) = (mabCosh6',mabsinh6') , (2) 

^As an example, for the Ising ferromagnet these specifications amonnt to consider a temperature slightly below 
the critical value Tc, in absence of external field. 
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Figure 1: A uniform boundary reflects a low energy kink. 

where niab is the kink mass (inversely proportional to the bulk correlation length) and 0 is 
known as rapidity. Two vacua |0a) and lO;,) (as well as the corresponding phases) are said to be 
adjacent if they can be connected by an elementary kink; when the connection requires a state 
IA'aci (^i)Arcic 2 (^ 2 ) • • • with n necessarily larger than one, the two vacua are said to 

be non-adjacent. 

As a further step towards the study of the wedge problem, we consider this statistical system 
on the half plane x ^ 0. We call boundary condition of type a a uniform (i.e. y-independent) 
boundary condition at x = 0 favoring phase a in the bull<§|, in such a way that the order 
parameter approaches {a)a as x ^ +oo. We will use the notation |Oa)o for the vacuum state of 
the quantum system on the half line with this boundary condition; more generally, the subscript 
0 will be used to indicate the presence of the vertical boundary. 

Phase separation can be induced through a change of boundary conditions. Within the held 
theoretical description, the change of boundary conditions from type a to type 6 at a point y on 
the boundary is realized by the insertion of a held fj,abi0,y), with non-zero matrix elements on 
states interpolating between |Oa)o and |Ob)o. When these two vacua are adjacent, which is the 
case we consider in this section, the simplest matrix element of yab 

o{OaMO,y)\Kba{0))o = e-™^“^"Vo(0), (3) 

where fo{9) gives the amplitude for the emission/absorption of a kink from the boundary con¬ 
dition changing point. The kink travels towards the boundary for 9 < 0 (in-state), and away 
from it for 0 > 0 (out-state). In- and out-states are related by the scattering operator [M] . 
As we are going to see, our computations involve low energy particles, whose scattering with 
the boundary is necessarily elastic, i.e. conserves the number of particles. Moreover, the held 
flab acts on a uniform vertical boundary, which preserves the energy. For these reasons the low 
energy scattering of a particle on the boundary is a pure rehection (Fig. [T]), and the relation 
between the in- and out-state for ([3|) takes for small momenta the simple form fo{9) = ±/o(—0). 
On the other hand, only the choice 

fo{9) = -fo{-9), 9^0 (4) 

^In a ferromagnet this is achieved applying a magnetic field on the boundary. 

®Here and below, in order to simplify the notation, we drop the indices on the kink mass. 
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Figure 2: Wedge geometry with boundary condition changing points at (0, ±i?/2) and an inter¬ 
face running between them. 


implies the property /o(0) = 0 which will eventually be responsible for the impenetrability of 
the wall. Generically, we will then have fo{0) = ci9 + 0{9^), with ci a constant. 

Passing from a vertical boundary to one forming an angle if: with the vertical involves a 
rotation in Euclidean space, and then a relativistic transformation for the associated quantum 
field theory. Recalling ([2]), this transformation shifts rapidities by so that the kink emission 
amplitude in the rotated frame, that we denote by /^, is related to that in the original frame as 

U{9) = fo{0 + i-ip); (5) 

our considerations on /o then yield 


/^(e) ~ ci(6»-hzV'), |(9|, |V^|<1. 


( 6 ) 


At this point we are able to consider, instead of the half plane, the more general wedge 
geometry of Fig. [2l The points (0, ±i2/2) are boundary condition changing points, such that 
phase b (resp. a) is favored for \y\ > R/2 (resp. \y\ < R/2) on the wedge. For mR large, i.e. 
when the system is observed on a scale much larger than the bulk fluctuations, one then expects 
an interface running between the points (0, ±i?/2), separating an inner phase a from an outer 
phase b. These expectations emerge from the theory in the following way. For \y\ < R/2 the 
order parameter in the wedge, that we denote by {a{x,y))wtab^ reads 


{a{x,y))w,^. 


(/)(06|/i.ba(0, 2 )^(^) y)A*ai)(0, 2)|06)— 


^Wbab 


(7) 


where the subscripts ±"0 indicate the different rotations performed for positive and negative y, 
and 


^Wbab = ip{0b\b‘-ba{0,R/‘2')yab{0,—R/2)\0b)--,p 

~ f + (8) 

This result is obtained inserting a complete set of particle states in between the two fields, taking 
the limit mR large which projects on the lightest (single-kink) intermediate state and to small 
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rapidities, and considering '0 small in order to use ([6]); here and in the following the symbol 
indicates omission of terms subleading for mR large. In a similar way we obtain 




^Wiab J-oo (27r)2 

X U{ei){KU0i)\<yM\Kab{e2))f-^{e2), (9) 


where we evaluate the order parameter field on bulk states, implying that the boundary condition 
changing fields account for the leading boundary effects at large R. The matrix element in Q 
contains a disconnected part proportional to 5{6i —O^) which yields a constant after integration, 
and then does not contribute to the derivative with respect to x we are going to take in a moment; 
the behavior of the connected part in the relevant region 61,62 0, is instead determined by 

the ‘kinematical’ pole (see |16] and references therein) 


{Ki,a{ 6 i)\a{^,{))\Kah{ 62 ))connecteA ^ * 


,(o-)b - 


61-62 


01 ~ 02 . 


With this information we obtain m 


dx{(T{x,y))wba 


{cr)b - {cr)c 


8 V 2 


I (x + ^) - 

RJ (1 + mRip"^) 


=-x 


where 


= VT^, 


e = 


2y 


X = 


2m X 


R ' ^ \l R K 

and, integrating back over x with the condition {a{+oo,y))wbab ~ (^)bi 


{(T{x,y))wba.b ~ (^)a + [{(^)b - (o-)a] 


erf(x)-^ 


2 X + y/ 2 mR ^ 


TT 1 + mRip"^ 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


for V' = y = 0 and {a)a = —{(^)b this result coincides with that obtained in m from the lattice 
solution of the Ising model on the half plane. 

It was shown in |13l [H] that the leading large R contribution to the order parameter profile, 
i.e. the one associated to the pole in (fT0]l . corresponds to a sharp phase separation between pure 

phases. In the present case of adjacent phases, there will be a single interface, with a probability 

pW 

(x;y) to intersect the line of constant ordinate y in the interval {x,x + dx). It follows that 
the leading large R expression of the order parameter can be written as 

p ^ p 00 

W{x,y))wiab ^ {^)b duP^'^\u;y) + {a)a du y), (14) 

J X J X 

where x{y) is the abscissa of the point on the wedge with ordinate y; this expression shows that 

p<P) 

(x;y) actually coincides with (fTTT) . Also for later use we introduce the additional notations 


A = 



X 




(15) 
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Figure 3: Contour plot of the passage probability density P^\x]y)/m for mR = 25, if: = 0.04. 
The leftmost contour line corresponds to P^\x',y) = 0, and then to the wedge. 

and rewrite this result as 

P?\x-y) ~ ^ + ■ (16) 

VTTK-^All + 

The requirement dxP^^\x‘,y) Ri 1 is satisfied as long as y/mR'ijj <C 1. Notice that (fTTl) or 
m show that P^\x;y) vanishes for \y\ = ^ + -fj which for the present case of small "0 are 
the coordinates of the wedge {x > —R'tp/2)-, hence, the properties (l5|), ([ 6 |) that we identified in 
momentum space indeed lead to an impenetrable wedge in coordinate space. A plot of (x; y) 
is shown in Fig. [3l 

3 Third phase and double interface 

In this section we still consider the wedge geometry of Fig. [2] with the same boundary conditions 
bab of the previous section, but now we study the case in which the phases a and b are not 
adjacent. More precisely, we consider the case in which the lightest state connecting |0a) and 
|0fe) is the two-kink state \Kac{0i)Kcb{ 02 )) , with a unique choice of the intermediate vacuum |0c). 
This situation arises, in particular, in the (dilute) ( 7 -state Potts model at hrst order transition 
that we discussed in the introduction. Indeed, the model is exactly solvable (integrable) in 
the scaling limit, and it is known that there are no kinks directly connecting two ferromagnetic 
vacua at the hrst order transition [351I29| . The lightest state connecting two such vacua |0a) and 
| 0 fe) is the two-kink state \KaoKQb) passing by the disordered vacuum |Oo); the symmetry under 
permutations of the q ferromagnetic phases which characterizes the Potts model m ensures 
that the elementary kinks Kao, K^a (o = Ij ■ ■ ■ j O') all have the same mass m. 

Technically, the difference with respect to the previous section is that now the large R 
expansion of d?]) is dominated by the contribution of the two-kink state. In particular, the 
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relevant matrix element for the boundary condition changing fields is no longer ([3]) but 

o(Oa|//a6(O,i/)|i^be(0l)i^ca(02))o = g—/ q(01, @ 2 ) ■ (17) 

As before, the low energy scattering properties of the kinks on a vertical wall can be used to infer 
properties of the amplitude /o(0i, 02 ), but now we will also exploit the integrability of the scaling 
(dilute) Potts model at the first order transition. Integrability ensures that the interaction of the 
two kinks on the wall can be regarded as consisting of two independent processes (factorization 
of the scattering [36]), and this in turn allows us to write a relation like (HI for each particle, i.e. 

/o( 01 , 02 ) = —/o( — 01 , 02 ) = —/o( 01 , — 02 ) , 01,02 —)-0. (18) 

Integrability also yields the exact bulk scattering matrix of the scaling (dilute) Potts model at 
the first order transition |35l I29| . From this one reads, in particular, that at low energy the 
state |A'bo(^i)7iloa(02)) scatters in the bulk into the state — 1 ^ 160 ( 02 ) 1 ^ 00 ( 01 )), so that we have 
the additional relation 

/o(01,02) = -/o(02,0l), 01 , 02 ^ 0 . (19) 

Equations (fTSl) and (fT^ lead to 

/o( 0 l, 02) ^C2 0102(0? - 0 i), 01,02 < 1 . ( 20 ) 

As before, the passage from the vertical boundary to that rotated by an angle involves a 
rapidity shift, 

^( 01 , 02 ) =/o(0i+#, 02 + *iA), (21) 

and the leading large mR expression for the order parameter in the wedge can be written a^ 

~ '^^^yy04 /v.(g4,03)(i^6c(03)l^ca(04)|a(O,O)|iPoc(0l)i^c6(02)) 

X /-v.(0i,02)m#2,03,04;x,y), (22) 

where 

3^(01,02, 03, 04; a;, y) = I7“(0i;x,y)C/“(02;x,i/)I7+(03;x,y)I7+(04;x,y), 

U^{9'X y) = e™(-f ±y) ™®h 0 =Fimxsinh 6 » 


(23) 

(24) 




f d0id02 

JrI (2vr)2 

3c? 

vr (mR)^ 


/v.(02,0i)/-^(0i,02)e-^(“^^^i+“^^''^) 




(25) 


The matrix element of the field a on two-kink states entering (I22p contains three types of 
contributions, depending on the number of annihilations that arise when particles on the left 
and on the right have the same rapidity. The three contributions are schematically depicted in 

'^Generically, we keep the notation c for the third phase; c = 0 for the Potts case. 
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Figure 4: The four-leg matrix element of the order parameter field a decomposes into the sum 
of the connected and disconnected parts in the r.h.s. 


Fig. [Hand correspond to a connected part (no annihilations), a partially disconnected part (one 
annihilation) and a totally disconnected part (two annihilations). The complete computation of 
the order parameter taking into account all these contributions has been performed in [T4] for 
the case of the strip geometry and in m for the case of the half plane {ip = 0). As expected 
from the fact that the two-kink state yields the leading contribution, the results correspond 
to the presence of two interfaces separating the intermediate phase c from the phases a and b. 
For the case of the wedge the complete calculation becomes cumbersome, and it is particularly 
interesting that we can still obtain the complete results in a relatively simple way through the 
following procedure, whose exactness we explicitly checked for the strip and the half plane. 

Generalizing what seen in the previous section, the probability P 2 ^\xi, X 2 ',y) that one in¬ 
terface intersects the line of constant ordinate y in the interval (xi, xi -|- dx), and that the other 
interface intersects the same line in the interval (x 2 ,X 2 + dx) is related to the order parameter 
as 


where 


{a{x,y))w,,k = / 

J. 

r+oo 

dxi 

i 

/■+°o . 

/ dX2P2 {xi,X2]y)(T{x\xi,X2) , 

J X 

(26) 


(o-)a, 

X < X < min(xi, X 2 ), 


cr(x|xi,X2) = < 

(o-)c, 

min(xi, X 2 ) < X < max(xi, X 2 ), 




X > max(xi, X 2 ), 



and x{y) is the abscissa of the wedge. On the other hand 


roo 

{xi;y)= / dx2 (xi, X2; y) 

J X 


(27) 


is the probability that one of the two interfaces passes in the interval (xi,xi -|- dx) at ordinate 
y, irrespectively of the other. Since it is the held a which ‘detects’ the interfaces, it is natural 
to expect, and we checked explicitly that this is the case for the strip and the half plane, that 
Pi^{x',y) is determined by the second term in the r.h.s. of Fig. |H proportional tcl§ 


{Kba{03)\(r{0,0)\Kab{ei))connectedS{92 - 9^) , 


(28) 


“Of course there are analogous terms with different pairings of rapidities, all giving the same contribution to 


(|22l). 









and to which we refer as the two-leg term, from the number of particles connected to the field 
a. Up to the factor <5(02 — ^ 4 )) corresponding to the undetected interface, this two-leg term is 
the same we studied in the previous section for the single interface. Plugging (1281) into (j22p we 
obtain 


{cr{x,y)) 


two-leg 


OC [ d0id02d03d04/^(04,03)^?^^/-^(01,02) 

jRi til — (73 

X g-™[f Et=l + f ^ 


(29) 


where, as usual, we took into account that small rapidities dominate at large R and we used 
(HOD; as in the previous section, a derivative with respect to x cancels the pole. On the other 
hand, we can also write ((t(x, in a way analogous to ([HD, with P ^2 replacing P^\ 
and ((t)c replacing {a)a or {(j)h- It follows that ^ i-®- 



OC [ dx2 [ d6'id6»2d6»3d6'4/.0(6(4,6»3)/_^(6'i,6»2) 

J Jr* 


(30) 


where we used 5{z) oc f dse*^^. Comparison with (|27l) shows that the integrand of the integral 
in X 2 in (l30l) is proportional to P^^^. Since the identity 

[ d/3id/32/v,(/3i,/32)e-^+'5i^i+''?2^2 =27r/_,^(gi,g2)e-^, (31) 

Jr^ 

holds for the function defined by (1201) and (1211) . we finally obtain 


P2^’^^(xi,X2;y) 


l^/_i(i+,)^(f/i,r? 2 )/_i(i_,)^(??i,r? 2 )e ^2 ( 32 ) 

N~ 

J 2 ^ + (1 + t)ii, r?2 + (1 + €)'ip) fo{rji + {l- €)ip, ?72 + (1 - e)^)e"^i"^2 , 


where we are using the notations (1121) and (jlSp with r]i and Xi corresponding to Xj, and is 
dimensionless and determined by the condition dxPi^ 2 ix] y) = 1. 

Recalling the form (I20p of the function /o, we see that the joint passage probability density 
(PI21) factors the terms {rji — 772 )^ and ±e = 1 -|- rn/ili (i.e. Py = ^ + ^)- It follows that the 
considerations in momentum space that led us to the result m for the function f^{ 01 , 62 ) 
produce in coordinate space a mutual repulsion among the interfaces {p!^\x,x] y) =0), as well 
as the presence of an impenetrable wedge along which the passage probability density vanishes. 
A plot of P ^2 {p y) is shown in Fig. [S] 

For V’ = 0 (I 32 D reduces to 


p!f\xi,X 2 \y) 


idxbdhdz 2 ^^-xi-xi 

37r k 2 A 2 


(33) 


a result which is known [38] to correspond to the so-called “vicious” walkers [33] on the half line 
X > 0: the walkers start at x = 0, move randomly with the constraint of avoiding each other 
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Figure 5: Contour plot of the passage probability density m = 25 and 

Ip = 0.04. The leftmost contour corresponds to = 0, and then to the wedge. 



Figure 6 ; Order parameter M“^(cj 3 (x, and contour lines for the dilute three-state Potts 

model at first order transition. Due to permutational symmetry, <73 does not distinguish between 
phases 1 and 2 , and the intermediate bubble of the disorderd phase is clearly visible. 
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Figure 7: Order parameter M~^ {a{x,y))^^ ^ for the Ising model in the undilute (continuous) 
and dilute (dashed) cases; in the dilute case the intermediate disordered bubble smoothens the 
profile. 

and the boundary, and return to x = 0 after a time R. Hence, our result (j,S2l) yields the exact 
joint passage probability density of two vicious walkers in the wedge. 

The order parameter can be determined from through (1261) . We quote here the explicit 
result in the case of the half plane, for which the expressions simplify. For the Potts model the 
order parameter held has components ak {Ylk=i = 0), with {ak)o = 0 in the disordered phase, 
and we obtaiij®! 

(crfc(x, = [{(^k)b + (o-fc)a] A{x) - 2{ak)aB{x) + (o-fc)a , (34) 

with 

(o-fc)a = a,/c = l,...g, (35) 

Mx) = --^X^ (x^ - 3 ) e“2^" + ^j^X (-2x^ + X^ - 6) e"^"erf(x) + exi{xf , (36) 

B{x) = (-6 + X^ - 2x'^) + erf(x). (37) 

dvvr 

A plot is shown in Fig. [6] for q = 3] for q = 2 (Ising) Fig. [7] compares the order parameter prohle 
in the dilute case with the undilute result of Ref. |14) . 

4 Conclusion 

In this paper we developed the theory of phase separation in a planar wedge for the case in which 
a macroscopic bubble of a third phase forms in between the two phases favored by the boundary 
conditions. We discussed the g'-state Potts model (dilute for g < 4) at its hrst order transition 
as an example to which the theory applies. In principle the full held theoretical calculation is 

®These results coincide with those derived by direct summation of all terms in Fig. |4][37]. 
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much more complicated than that performed in m for the case of a single interface. However, 
we found that, isolating a specific contribution to the order parameter which corresponds to the 
detection of a single interface, the formalism allows to reconstruct the complete result. This 
finding, that we checked explicitly against the full calculations for the cases of the strip and of 
the half plane, appears promising for further developments. 

For the case of the wedge, it is worth stressing that the very fact that the final result exhibits 
a wedge-shaped path along which the passage probability for the interfaces vanishes provides a 
non-trivial consistency check for the theory. Indeed, the calculation starts from considerations 
in momentum space, in which the presence of a boundary is codified in properties of matrix 
elements of boundary condition changing fields. Moreover, the boundary is initially flat, and 
the information about the wedge is introduced through relativistic transformations performed 
on the matrix elements, always in momentum space. While in principle these are the same logic 
steps performed in m for the separation of two phases, in practice the present case with a third 
phase is much more structured and leads to the specific form (I2UI) . (1211) for the matrix element 
of the boundary condition changing field. It is then remarkable to realize going on 
with the computation that the wedge in real space emerges because 6 * 2 ) turns out to fulfill 

the self-Fourier transform property (1311) . In this way, the mechanism which eventually accounts 
for wedge covariance acquires surprising mathematical implications in presence of a third phase. 

This appears to have also additional implications. Indeed, we arrived at (l20|) exploiting also 
the integrability of the scaling Potts model, which in turn ensures factorization of the scattering 
and equation (fTSl) . On the other hand, since (fTSl) is necessary to arrive at (l20|) . and the latter is 
necessary for the appearance of the wedge in real space through (I^TT) . we are led to conclude that 
factorization of the scattering at low energies is required in systems allowing for the appearance 
of a bubble of a third phase. Notice that this is a weaker property than integrability, which 
implies factorization of the scattering at all energies. 

We determined the joint passage probability for the interfaces separating the three phases 
and found that, in the case of a flat boundary (tilt angle 'ijj = 0 for the wedge), it coincides with 
the known probability for vicious walkers in the half plane. Hence for ip ^ 0 our result provides 
the passage probability for vicious walkers in a wedge. The name vicious walkers is used in the 
literature for random walkers subject to the constraint of avoiding each other and the boundary. 
In our framework the properties of random propagation and avoidance for the interfaces emerge 
from the limit of large separation R between the boundary condition changing points on the 
boundary (pinning points for the interfaces), which is needed to observe phase separation. This 
limit projects the dynamics of the particles to low energies, where it turns out to reduce to 
fermionic statistics and becomes universal. Consistently, the average distances between the 
interfaces and between the interfaces and the boundary grow as and are much larger 

than the range 1 /m of the particle-particle and particle-boundary interactions, whose details 
then affect only subleading orders in the large mR expansion. 
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